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Abstract 

The known Peres-Horodecki criterion and scaling criterion of separability are considered 
on examples of three-mode and four-mode Gaussian states of electromagnetic field. It is 
shown that the principal minors of the photon quadrature dispersion matrix are sensitive to 
the change of scaling parameters. An empirical observation has shown that the bigger the 
modulus of negative principal minors, the more entangled the state. 
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1. Introduction 

The entanglement phenomenon is one of the most important quantum properties of multimode 
Gaussian states (both pure and mixed) was studied in the context of their separability and 
entanglement in [U [2] where the Peres-Horodecki criterion of separability [3], H] was implemented 
for detecting the entanglement of two-mode-light Gaussian state. The Peres-Horodecki criterion 
of separability of a quantum state is based on the property of positive but not completely positive 
map of the density matrix [5J, which is the transpose transform of the density matrix. For a 
separable state of a bipartite system, the partial transpose transform applied to variables of one 
subsystem yields a new density matrix that corresponds to the physical state. For the density 
matrix of an entangled state, the partial transpose transform can give a matrix that is a negative 
Hermitian matrix. The criterion (called ppt-criterion) was used in [H |2] to study the entanglement 
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of Gaussian states of photons. Later on, other researchers successfully applied the ppt-criterion 
(see, for example, [6-9]). 

The other separability criterion based on a nonpositive map of the density matrix, called the 
partial scaling transform, was suggested in [THl [EE] ; its implementation is connected with proper- 
ties of the Wigner function [12J with respect to the scaling of momentum. For a separable state 
of a two-mode light, the scaling of the second component of momentum provides an admissible 
Wigner function. The scaling criterion can be extended to the case of multipartite quantum states. 
The ppt-criterion is a special case of the scaling criterion with the choice of the scaling parameter 
equal to —1. The scaling criterion was used recently by Chirkin and Saigin [T31 [33] to study the 
entanglement properties of three-mode light generated in nonlinear crystals due to the interaction, 
where three different processes are involved, namely, two parametric down-conversions and one 
parametric up-conversion. 

In this work, we study the application of the scaling criterion to three- and four-mode photon 
states. We consider both the pure and mixed states. Our aim is to find some empirical correlations 
of the degree of entanglement and negative values of principal minors of the photon quadrature 
dispersion matrix as functions of the scaling parameters. 

The paper is organized as follows. 

In Sec. 2 the definition of separability is discussed. In Sec. 3 the Peres-Horodecki criterion 
is reviewed. The properties of scaling transform are discussed in Sec. 4. The application of the 
scaling criterion to the three-mode case is elaborated in Sec. 5, while the four-mode case is studied 
for a pure Gaussian photon state in Sec. 6. A mixed four-mode Gaussian state is considered in 
Sec. 7. Some conclusions and perspectives are presented in Sec. 8. 



2. Separability and Entanglement 

We consider a single-mode photon state described by the quadrature operators q and p. The 
density matrix p of any field state obeys the following conditions: 

p + = p, Trp = l, p>0. (1) 

When taking a transpose of the density matrix, these conditions are kept for a new 'transposed' 
state: 

(P T ) + = P T , Trp T = l, p T >0. (2) 

The transpose transform of the density matrix is equivalent to the time inversion or complex 
conjugation of the density matrix. There is a one-to-one correspondence between density operators 
and Wigner quasidistribution functions W(q,p). From the definition of the Wigner function it 
follows that the transpose transform of p is equivalent to a mirror reflection of momentum p in 
the phase space, namely, 

P — P T W(q,p) W(q, -p). (3) 
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Now we consider a two-mode state with the density operator 2) and annihilation operators 

„ _ qi+ipi „ _ g 2 + ip2 

We say that a state is simply separable if 

p(l,2)=p(l)®p(2). (4) 

A separable two-mode state is defined as a state with the density operator representable as a 
convex sum of simply separable states: 

p(l,2)=X;*i# ) (l)®P? ) (2), £% = 1, fci>0. (5) 

3 j 

Here pi^(l) and ffi\l) are the density operators of modes 1 and 2, respectively. The states which 
cannot be presented in such a form are called entangled states. In other words, by definition, an 
entangled state is a state which cannot be presented as a convex sum of simply separable states. 



3. The Peres- Hrodecki Criterion 

It is evident from (4) that the partial transpose, i.e., the transpose of a matrix for the second-mode 
term pi\2) — > (p2^(2)) , leads to a separable positive density operator 

^(1,2)=E^P? ) (1)®(#(2)) T (6) 
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Thus, p p *(l,2) describes a new separable state. 

Summarizing, we can say that the partial transpose p — > fp 1 of a separable state necessarily 
gives the density operator, which fits the conditions (1). This is the Peres- Horodecki criterion of 
separability. 

Rewriting the decomposition (5) in terms of Wigner functions we obtain 

W( qi , Pl ,q 2 ,p 2 ) =J2 l 3W?\qi,Pi)W?\q 2 ,p 2 ). (7) 
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Let us consider the Wigner function of the generic Gaussian form 

W(q,p) = -±=exp (-Iq^Qt) , (8) 
where the four-dimensional vector Q reads 

Q = (<?i - (gi> ,92 - (q2) ,pi - (pi) ,P2 - W) ( 9 ) 
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and the matrix a is a 4 x 4 real symmetric variance matrix 

1 



Or 



(10) 



with fi = gi, f 2 = g 2 , r 3 = pi, f 4 = p 2 . 

From (3) and (8) it is easy to see that the partial transpose is equivalent to the transform of 
a variance matrix a — ► a', which is defined as follows: 
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The conditions of existence of a partial transposed state in the Gaussian case can be rewritten in 
the form of the Robertson-Schrodinger uncertainty relation [T5l HE]. Thus, the Peres-Horodecki 
criterion of separability is equivalent to the condition of positivity of principal minors of the 
matrix 
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(12) 



(13) 



Here the dimensionless units with H = 1 are used. 

Since the transform considered does not change the minors but the determinant, the Peres- 
Horodecki criterion can be simplified into the inequality 



det S > 0. 



(14) 



This criterion is a necessary and sufficient condition of separability of Gaussian states in the 
two-mode case, but the application of the Peres-Horodecki criterion to the n-mode case with 
n > 3 fails to detect all entangled states. 



4. Scaling Transform 

Now we consider a scaling transform, which is applicable for detecting the entanglement of Gaus- 
sian states. We generalize the reflection transform p — ► —p in the phase space and use the scaling 
of momentum p — > Xp. In that way, we define the map 

j p — ► Xp. 



4 



For A G [— 1, 1], this map defines a semigroup of maps. The map (15) is equivalent to the scaling 
of time t — > At In the case A = 1, this is an identical map; for A = —1 we have a partial transpose 
used in the Peres-Horodecki criterion of separability, which was discussed in the previous section. 

Upon applying the scaling transform, the single-mode Wigner function W(q,p) modifies as 
follows: 

W(q,p) — W x (q,p) = NW(q,Xp), (16) 

where N is a normalization constant. 

The initial Wigner function is normalized 

dqdp 



J W(q,p)'- 
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and the transformed Wigner function must be normalized, as well: 



dqdp 



! 



l. 



Substituting p' = \p into the integral (18), we obtain 

Taking into account (16) we have the following result for the normalization constant iV: 

N 

Hence, the scaling transform modifies the Wigner function as follows: 

W(q,p)^W x (q,p) = \\\W(q,\p). 



N = |A|. 



(17) 



(18) 



(19) 



(20) 



(21) 



Now we find out how the partial scaling transform affects the dispersion matrix of the two- 
mode quantum state. We perform the scaling of the second momentum by the factor of A, 
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(23) 



X"P2Ql X°V111 X a P2Pl X 2<J P2P2 / 

We obtain a test of separability based on the Robertson-Schrodinger uncertainty relation for 
the partially scaled matrix o" A 

det (a x + ^ > 0. (24) 
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5. Scaling Transform in the Three-Mode Case 



Changing the sign for one momentum only may not be efficient in the case of three and more 
modes. The Peres-Horodecki criterion detects the bipartite entanglement. This limitation can be 
overcome using the scaling transform by changing momenta for several modes in different ways. 
In the three-mode case, we can apply the following transform: 



(25) 



The scaling parameters are determined, as in the previous case, in the domain from —1 to 1, 
namely, Ai G [—1,1], A 2 G [—1,1], and A 3 G [—1,1]- The applied transform leads to the scaled 
variance matrix 
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The uncertainty relation for the scaled matrix demonstrates the criterion of separability. Since 
the first three principal minors remain unchanged, taking into consideration the initial uncertainty 
relation, the criterion reduces to the three inequalities: 



and 
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OAiAaAa + ^ 



C^Ai A2A3 + 



-/ 
/ 

'0 -/ 
/ 



CAiA 2 A 3 





I 







>o, 
>o, 

> 0. 



(27) 
(28) 

(29) 



Here / is the 3x3 identity matrix. The fifth and fourth principal minors are denoted as M 5 and 
M 4 , respectively. The criterion consists in checking the sign of functions S(Ai, A 2 , A3), S 5 (Ai, A 2 ), 
and S 4 (Ai) in the domain Ai x A 2 x A 3 = [—1, 1] x [—1, 1] x [—1, 1]. The negativity of one function 
is an indicator of entanglement of the given state with a certain variance matrix. 
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Below it will be shown that the absolute value of negative function £(Ai, A 2 , A3) is connected 
with the degree of entanglement. 

We consider a specific pure three-mode Gaussian state with the wave function as follows: 



ty(x,y, z) = N exp 



x 
"2 



y_ 
2 



+ cixy + c 2 xz + c 3 yz 



(30) 



where ci, c 2 , and C3 are the parameters of quadratic form determining the degree of entanglement. 
The normalization constant N reads 



N = 



c| — C3 — 2C1C2C3 
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(31) 



In this case, the variance matrix a TiT . can be calculated using the definition (10). After calculating 
the elements of a, we obtain the following result: 



a = 
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(32) 



where we use the notation 
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(33) 



V = 
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The uncertainty relations for the state under study and the normalization conditions determine 
the admissible domain of the coefficients of the quadratic form. The following relations show the 
possible domain of the coefficients: 



r cie(-M), 

c 2 G(-l,l), 

c 3 e(-l,i), 

1 - c? - c? - 



(34) 



Co — 2cic 2 c 3 > 0. 



The Peres-Horodecki criterion for the three-mode case is obtained by applying the transform 
of quadratures as follows: 

' Pi — >Pi, 

P2 — ► -P2, (35) 

I P3 >■ -P3- 

The Robertson-Schrodinger uncertainty relation for the transformed matrix provides two inequal- 
ities - one from determinants of the whole variance matrix and the other from the fifth principal 
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Fig. 1. Contour plots of £(Ai, A2, A 3 ) for the wave function (36), where Ai = 1/2 and c\ 
the left) and c\ = 5/6 (on the right). 



2/3 (on 



minor of the matrix. The determinant [analogous to (24) in the three-mode case] is equal to zero 
for any Ci, C2, and C3. The determinant of the fifth principal minor is proportional to cf, which 
leads to the conclusion that the Peres-Horodecki criterion does not detect any entanglement for 
all wave functions of the form (30) with c\ = 0. 

Below we examine the application of the partial scaling criterion to the wave function (30) for 
certain coefficients Ci, c 2 , and c 3 . 

First, we consider the case where the entanglement is given only by the term xy or, in the 
language of coefficients, the case 02 = 03 = 0. Obviously, the problem is symmetric with respect 
to permutations of ci, ci and C3. From relations (30) and (31) and conditions (34), we have the 
following form of wave function for any c\ G (—1,1) 
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Simplifying the variance matrix (32) and applying the partial scaling transform (25) we obtain 

1 c?(l-A 1 A 2 ) 2 -(l-A?)(l-Ai) 
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(37) 



64A?A|A| 

For Ci = 0, expression (37) is nonnegative in the domain Ai x A 2 x A3 = [—1,1] x [—1,1] x 
[—1,1]. So, for this separable state, the uncertainty relation holds and the criterion does not 
show any entanglement. Otherwise, for any c\ G (—1, 1), such a pair Ai and A2 exists for which 
S(Ai, A 2 , A 3 ) is negative. Thus, the partial scaling criterion (27) detects the entanglement in this 
case. Moreover, it is easy to see that the rate and area of negativity increases with the coefficient 
Ci. These facts will be demonstrated by some examples below. 
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Fig. 2. Contour plots of S(Ai, A2, A3) for c\ = c 2 = 1/4, where Ai = 1/2 and C3 = 1/4 (on the left) 
and C3 = 1/2 (on the right). 



An illustrative example is shown in Fig. 1. 
Contour plots present the function £(Ai,A 2 ,A 3 ) 
with Ai = 1/2. To show that the graph is mov- 
ing down with increase in entanglement, we con- 
structed contour plots for Ci = 2/3 and 5/6. The 
negative area of the graph for bigger coefficient C\ 
is larger and lower than for smaller c\. Since we are 
not interested in the positive part of the graphs, the 
grey scale shows only the rate of negativity. The 
negative part of the graphs is bounded by lines, 
where the function is equal to zero. The area be- 
tween lines is a positive part of the function. It 
is worth mentioning that for c\ — *■ 1 or, in other 
words, for the maximum entanglement in this case, 
the positive part collapses to the line A2 = Ai. In 
the considered case, it is the line A 2 = 1/2. 

The completely entangled state leads to the 
same result. We demonstrate it by the case c\ = c 2 = 1/4. In Fig. 2 graphs for c 3 = 1/4 
and 1/2 are shown. Obviously, the graph for bigger C3 has larger negative area and its absolute 
value in points, where both functions are negative, is bigger. The zero level is shown by a black 
curve. The area inside it contains only positive values. 

A vivid illustration that the graph is mowing down with increase in entanglement is shown in 
Fig. 3 where the function is constructed for the wave function with c 2 = c 3 = 1/4. As before, it 




Fig. 3. Contour plot of S(Ai,A 2 , A3) for 
c 2 = c 3 = 1/4, where Ai = 1/2, A 2 = 1/4. 



represents E(Ai,A2,As) with Ai = 1/2 and A2 = 1/4. Obviously, the positive part of plot (white 
color) diminishes with increase in c\ from to 7/8 and the negative part is moving down. 



6. Scaling Transform in the Four-Mode Case 



To apply the graphing method discussed in the previous sections to the four-mode state, we apply 
the transform of momenta of the form 



(38) 
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Fig. 4. Contour plots of £(Ai, A2, A3, A4) for the wave function (40), where Ai 
and c\ = 1/8 (on the left) and c\ = 1/2 (on the right). 



A 3 

-1 and A 2 = 1/2 



All scaling parameters belong to the interval [—1,1]. The transformed variance 8x8 matrix 
0"a 1 a 2 a 3 a 4 is similar to (26). The Robertson-Schrodinger uncertainty relation can be easily reduced 
to four inequalities, in complete analogy to the previous case. We consider only the inequality for 
determinant of this matrix. It will be shown that the function defined below gives results very 
similar to the ones obtained in the three-mode case. The inequality under consideration reads 



S(Ai, A 2 , A 3 , A 4 ) = det 



OA1A2A3A4 + 




I 







> 0. 



(39) 



It is worth mentioning that the Peres-Horodecki criterion applied to the four-mode wave 
function does not give complete information on entanglement - in some cases, it does not work. 
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We investigate wave functions of the form analogous to the ones considered in the three-mode 
case. First, let us look at the function 

*(x, y, z, k) = iVexp (-^ - V - - Z - - k — + Cl xyj , (40) 

where 



'1 — <r 

N = 1 1 . (41) 

7T 

Graphs of E(— 1, 1/2, A3, A4) for c\ — 1/8 and 1/2 are shown in Fig. 4. One can observe the 
same effect of moving the graph down with increase in entanglement. Entire graphs are located 
in the nonpositive half-space. 

For completeness, in Fig. 5 contour plots for completely entangled states are shown. The 
graphs are constructed for 

, f x 2 y 2 z 2 k 2 1 1 11 1 , \ , . 

y, z, k) = Nexp I -— - — - — - — + -xy + -xz + -xk + -yz + -yk + c 6 zk I , (42) 

where c 6 = 1/8 and 1/2. 

In both cases, the partial scaling transform detects the entanglement. Also it can be seen 
that the bigger the entanglement, the lower the location of the plot of £(Ai, A2, A3, A4). The two 
functions are nonpositive in the domain A3 x A4 = [—1, 1] x [—1, 1] too. 
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7. Scaling Transform Criterion for Mixed Gaussian States 



In this section, the partial scaling method is tested on the applicability for mixed Gaussian states. 
First, we consider the three-mode state with the variance matrix 
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(43) 
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The graph for S(l, A2, A3) is presented in Fig. 6 (on the left) where the positive part of the 
graph is shown with white color. Obviously, the partial scaling method detects the entanglement 
in this case - the graph has a negative area in the domain A2 x A3 = [—1, 1] x [—1, 1]. 





A, 



Fig. 6. Contour plot of S(Ai, A2, A3) for the variance matrix (43), where Ai = 1 (on the left) and 
contour plot of S(Ai, A2, A3, A4) for the variance matrix (44), where Ai = —1 and A2 = 1/2 (on the 
right). 

A similar result is obtained for the four-mode mixed Gaussian state. In Fig. 6 (on the right) 
the graph for E(— 1, 1/2, A3, A4) is shown. 
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The 8x8 variance matrix reads 
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(44) 



In both plots, there is a negative area, which depicts the entanglement of states with the 
variance matrices (43) and (44). 



8. Conclusions 

To conclude, we point out the main results of our work. 

We have shown that the new criterion of entanglement based on the scaling transform of 
the photon quadrature components in the Wigner function of the quantum state provides the 
possibility to detect the entanglement for three- and four-mode fields. The calculated principal 
minors of the scaled quadrature dispersion matrices shifted by the quadrature commutator matrix 
become negative for entangled Gaussian states. The values of modulus of the negative principal 
minors can serve as empirical characteristics of the degree of entanglement. The results obtained 
for the three-mode entangled Gaussian states, in general, agree with the previous results of [131 [H] 
where the scaling criterion of entanglement was applied to the photon states created due to the 
light interaction with nonlinear crystal. In a future publication, we will extend the analysis of the 
entanglement based on the scaling criterion to the multimode light beams (both of Gaussian and 
non-Gaussian tipes) and compare this criterion with other approaches to the problem of detecting 
the entanglement. 
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